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MECHANICS BY QUATERNIONS. 



BY PROF. B. W. HYDE, UNIVERSITY OP CINCINNATI. 

(21). Center of Gravity. — We have seen (Eq. 21) that the vector of the 
center of mass of a system of particles whose weights are w lf w 2 , w s , etc., 
is given by the equation 

o — I( - W '^ 

' p0 — ^f 

To apply this formula to finding the center of gravity of any line, sur- 
face, or volume we have only to substitute for w the element of the line, 
surface, or volume multiplied by a coefficient of heaviness, aud then to sub- 
stitute integration for summation. Let h be the weight of a unit of volume 
of the substance. If it is the same at all points of the body it may be taken 
outside of the sign of integration, and will then cancel out ; if it is not con- 
stant it must be expressed as a function of t p before the integration can be 
performed. 

(22). Center of Gravity of a Line. — By this is meant the c. g. of such a 
body as a fine wire bent into some curve, the diameter being so small com- 
pared with the length that it may be regarded for practical purposes as a 
true mathematical line, i. e., as if the wire, still retaining its weight, should 
have its diameter reduced to nothing. Let the equation of the curve be 

ip = i? (*)• • • • d ip = y (o*> 

and ds = Td t p = T^'fy.dt — element of curve. Therefore, for this case, 

D _fphTd /P _ fh / <p(WW(t).dt .„. 

'^ fhTd lP fhT t (f'(t).dt { > 

(23) Center of Gravity of a Plane Area. — Let the bounding curve be 

iP = Mfi> 

and let e be some unit vector in the plane of the curve. We may take either 
1 2 3 4 as the element of area or 1 2 5 6. In the! 
first case we have area 12 3 4= TVs J pSe~ 1 d / p\ 
(at the limit) ; so that the area between a pair of | 
perpendiculars to s will be 

fTVe tP Se-id t p = j TVstfSe-* rf.dt, 

where for brevity t <p is written for t <p(t). Now the vector to the c. g. of the 
element is ,p + ^Vs t p = ,<p -\- ^eVetf; hence for the center of gravity of 
the area between e and the curve 

/Po = jh{j<p+le Vs / f)TVe,<pSe- 1 ,<p'.dt-i- fhTVetfSe- 1 t f'.dt. (57) 
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In order to get the element 1 2 5 6 we must evidently change the tensor 
of t p independently of the versor. Therefore write the equation t p = wtfif). 
By giving u a succession of values we shall obtain a series of curves such as 
5 6. If we differentiate with respect to u only T t p varies, so that the end 
of t p moves say from 1 to 6, while if we differentiate with respect to t the 
end of ,p moves along such a curve as 1 2 or 6 5. Let 

d,P — j) and d ' p — D 

then the area of element 12 5 6 will be TVD m D t .dudt. But X>„ = ,<p, and 
B t = utf', by the equation assumed above. Therefore 

lPt = ffhtfTVtftf'.tfdtdu -+-ffhTV,<p,<p'.udtdu; (58) 

or integrating for u from to 1, to cover the space from the origin to the 
bounding curve, 

//>. = hfhtfTVw'Jt ■+- ifhTV^tp'.dt. (59) 

This integration however can only be performed when h is not a function 
of u. 

Equations (58) and (59) hold equally well when ,p = t f{() does not rep- 
resent a plane curve. They give in any case the c. g. of the surface swept 
over by the radius vector. 

(24). Any Surface. — Let the the equation of the surface be l p = l f(x,y). 

Then D, = -£-, and D v = -*£- are vectors II to tangents to the surface 
dx dy 

at the end of ,p, for by supposing x to be constant we have one curve on the 
surface, and by supposing y to be constant we have another. By giving 
successive constant values to x and y the surface will be divided up into 
four-sided figures which when small enough may be regarded as parallelo- 
grams. The area of such an elem. parallelogram will be TVD x D v .dxdy. .• . 

/Po = ffh /P TVI> x D y .dxdy-*- ffhTVD x D y .dxdy. (60) 

(25) Any Solid. — If the equation of the bounding surface is ,p= t <p{x,y), 
then to obtain a parallelopipedical element we must vary the tensor of ,p. 
Therefore write, as in Art. 23, ,p = u t <p{x, y). By giving u successive val- 
ues differing by du, and extending from to 1 we divide up the solid into 
a series of shells. The solid element will then be SD K D x D y .dudxdy. But 
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If h is not a function of u we can integrate for u between and 1, thus 
obtaining 

A= i/M ^.dvdy+lffhS^^.dvdy. (62) 

(26). Surface of Revolution. — We will now apply eq. (60) to the general 
equation of a surface of revolution. Let ,p = ,<p(t) be the equation of any 
curve, plane or tortuous. If this curve be revolved about an axis along 
which e is a unit vector, and the origin be taken at a point of this axis, a 
surface of revolution will be generated whose equation will be 

1 — 1* 

This will appear from the fact that if be constant in eq. (63), the eq'n 
represents the generating curve turned through the angle 6 about e, while 
if t be constant the equation is that of a circle generated by the extremity 
of ,<p(t) revolving about e. 

To apply (60) we have to evaluate TVD x D y . 

20 e _e 

B x = Be = e "" Ve,f and D y = D t = t* ,<p't n , 
in which the t is omitted for brevity. 

.-. TVDJ) y = TVDeDt = 2T.t2*+* Fe,p.t*+*y «-*+■« 

= ZT[«7e,?&,p'oos0— Ve / fSe / <p'smd+eS.e / <p'Ve,f'] 

= -i/(S 2 .£^' Vs,<p— V%<pS%f)= v /( /<P ' 2 VH^—Sh^,?') 
= TV. t <p' Vttf. Therefore (60) becomes 

lPo = J" fh t** # t ~***T V.y Vt,<p dt dd -i-ffhTV.tf' Ve,<p dt dd. (64) 
If h be not a function of we may integrate for 6 from to lf obtaining 

/j0o = fhW^Serf— (e 2e,H " r - l) Vt#\TV.rf Vetfdt+dfhTV.y VttpdL 

(65) 
If the integration be for a complete revolution so that 0, = 2jt, then the 
equation reduces to 

,p 9 = fhe-iSe^TV.y V S/ <p dt-t-fhTV.y Vt,ip dt. (66) 

If the generating curve lie in a plane passing through the axis we shall 
have Se,<p,<p' = 0, so that in this case TV. t (p' Vstf — T t <p' Ve t <p. 

(27). Solid of Revolution. — For the solid generated by the radius vector, 
when its extremity generates the surface of eq. (63), we have as in Art. 25, 

t p = us®-*-* tftys— 0-*-* . 

*See Tait's Quaternions, Art. 351 
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We need these to evaluate 

1 —1 %L it 
8D u DeD t = SDuVDeDt = u 2 Se* t <pe "T.**" Fe^e^e n , 

or by the value of VDeDt of the last Art. 

SDuBeDt = v?S(e- i Se,<p+s- 1 Ve J <pcoad+ Fe^sintf) 

= u?[S.e,<p' Vs t (p Sejf — V 2 e l <pSe,<p'~\ = — 8-tf/f' Ve,<p.u 2 . 
Hence substituting this value and that of t p above in eq. (61) we have 

"°° =fffhu s £",</>e~*S;<P / f Ve,<p.dtdOdu+fffhu 2 S. / <f> / <p' Ve y <p.dtdddu.{67) 

[Here, and through the remainder of this paper, for want of sorts, <f> is 
substituted for <p and for subscripts and indices th is written for 6."} 

If h be not a function of we can integrate as in last article from to 
1} obtaining 

l p = ffhu^Vd^Se^— in?— IjVs^S.tfyVetf.dtdu 

-¥■6 i / ' fhu i S.,<p / f VtrfM du. (68) 
If 0j = 2n, for a complete revolution, 

,/),= e-ij'fhu 3 Se l <pS. l <p J fVs / </>.dtditr±-ffhii 2 S. / <f> / fVe l 4>.dtdu. (69) 

If the generating curve be in a plane through the axis, so that V. V t <j) ,<[>'- 
Verf = 0, we may replace the quantity 8.$$' Verf in eqs. (67), (68) and 
(69) by the equivalent expression TVtf^' Vs t <p. For with the above con- 
dition we have 

V 2 l <f> l fV%<p= PVtfrf'Verf = SKrftf'Vetf; .-. TVtftf'Vetf 

=*±S.tf t fV,tf. 

(28). Guldinus' Properties. — We have by Art. 26, the area of a surface 
of revolution formed by revolving a plane curve about an axis in its plane 
through an angle 2^-r-n, 

S = ^JT.fVerf.dt: 
also by eq. (56) omitting the h, i. e., supposing the density uniform, 

/Po =f l m'.dt+fT,f.dt. 

Multiply both sides of this eqn. by (27r^-n)/T / <p'.dt and operate by TVs; 
... ^TVe /Po .fT,fdt = ^fT,fVe,<f>.dt = S. (70) 

The factor outside the sign of integration in the first member is the nth 
part of the circumf of the circle in which the c. g. of the curve t p — ,<p{t) 
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moves, while the integral in that member is the length of the same curve; 
hence the area of the surface generated by a plane curve revolving about an 
axis in its own plane is equal to the length of the curve multiplied by the 
length of the path of the center of gravity of the curve. 

Again, the volume generated by revolving the area bounded by a plane 
curve about an axis in its own plane through an angle 2;r-f-w is, by Art. 27, 

and by eq. (58), omitting h as before, the c. g. of a plane area is given by 

/Po = ^j'f / <pTV,<p,fM t dudt+~ffTVrf / <p'.ududt. 
or clearing and operating by TV.e 

— TVe / p l) .ffTV / <prf'.ududt = —ffTV / <p / fVerf.u 2 dudt (71) 

= v. 

Hence the volume of the solid formed by revolving a plane area about 
an axis in its own plane is equal to the plane area multiplied by the length 
of the path of the center of gravity. 

(29) These properties may be extended to the cases of a plane curve or 
a plane area moving so that the center of gravity of the curve or area gen- 
erates some curve ,p= / <p(t) to which the plane of the curve or area is always 
normal. Let I be the length of the moving curve and A its area, S and V 
respectively the surface and volume generated and t p = ^(t) the path of 
the center of gravity. Then for the length of an elementary arc d t p the 
curve will coincide with the arc of its osculating circle for that point, and 
therefore Guldinus' properties will hold; hence dS = ITd^^lTrf'.dt, and 
dV = ATj<p'.dt; and by integration 

S -= if^f.dt; V = A.fT t f.dL (72) 

(30). A few examples to illustrate the application of these formulae will 
now be given 

Center of Gravity of the arc of a Cycloid. — If d be the angle turned thro' 
by the rolling circle, and a the radius of this circle; then the equation of 
the cycloid with the origin at a cusp may be written 

jP = ,<p(0) = ai(d — sind)-\-aj (1—cosO), 
in which i is a unit vector along the base, and j a unit vector along the tan- 
gent at the origin. Therefore djp-hdd = ^'(d) = ai(l — cos#) -f aj sin 6, 
and TjP(0) = a\/[2(l — cos/?)] = 2a sin $6. Therefore, by eq. (56), put- 
ting 6 for t, and supposing h constant, 

W°o = aj [*(0— sin0) +j{l— cos 0)] sin ±0.dd -=-J*sini0.d0. 
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By integrating this equation from to it, i. e. from the cusp to the ver- 
tex of the curve, we have ,p = %a(i-±j). If we go from the origin to the 
next cusp, the limits will be and lit, and t p = a(jra-j-|j). 
As another example take the tortuous curve whose equation is 
lP = a{it+W2.jf +$**>), 
and which is projected on the three reference planes into a common, cubic, 
and semi-cubic parabola respectively. d,p-^dt = j</>'(t) = a{i+jt\/2-\-M) ; 

.:T,f(t) = a(l+f). 

■•■ t p*= a$(it+lJt 2 ^2 + lW)(l + ?)dt-s- ${\+P)dt. 

If the integration be from to 1 we shall have 

For an example of the application of eq. (57) take the parabola whose 
equation is ,p = t at+e\t 2 , in which e is a unit vector along the axis, and p 
is a vector along the tangent at the vertex, so that Se,a = 0. 

.•. ^'(t) = ,a + et, Ve,<p = tVep ■— U,a, eVe t tp = — t,a, TVe t <p — at 
and Se* 1 ^' = t. Suppose h to be constant, then 

,fl = Cl(/a+eP)aP.dt-*- CaJPdt = &P£M 
J o" J o ft 

Let us apply eq. (59) to the equation of the parabola just used, now how- 
ever regarding ,a as a vector along any tangent and £ as a unit vector along 
the diameter through the point of contact of the same. Then 

TVrf,f = tTV(,a + l«Q(,o + et) = t 2 TV{ t ae — l,ae) = it 2 TV t ae. 

••■ w°o = f f\ t a+%et)?dt+ Cm = hat+iei?- 

This gives the c. g. of any segment of a parabola. 

The equation of a circle may be written t p = t (f>{0) = a(i cos d-\-jsm 6); 
.•. $> = a( — isin0+j*eos0), and TV t (p t <p' — a?. Therefore, by eq. (58), 
if h = constant, and 6 = t, 

y /) = a| I (icosd-{-jsuid)v?dOdu-+- I I udddu 
J —ay « J —ay « 

= \a{\— u»)[tsintf— joosflJtS-f- |(1— « 2 )[0]±£ 

2 1 — w 3 isinfl 

This gives the c. g. of a portion of a ring contained bet. two concentric cir. 
1— m 3 _ (1— u)(-l-\-u+n 2 ) __ 1+w+m 2 . . _2 1+w-H^ sing ■ 

i__ M 2 (i_ M )( 1+M ) - 1+tt ; ■ • /j°o 3 • 1+u d m. 

If u = 1, we have for an arc of a circle, /( o = (sin#-j-#).ia. 
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Ellipsoid. — The equation of this surface may be written 

in which a, b. c are the semi-axes. The c. g. of the solid ellipsoidal shell 
will be found first by eq. (61), and then that of the surface, from this by 
making the thickness of the shell infinitesimal, for the reason that the sca- 
lar in eq, (61) is much simpler than the tensor in eq. (60) in the case of the 
ellipsoid. We have then 

-£— = — fi sin x-\-fi cos x cos y-\-fi cos x sin y, 

-j£_ = — fi sin x sin y-\-fi sin x cos y. Whence 

S,d> -$--¥- = 8,a,b.c sin x = — abc sin x. 
ax ay 

Thus eq. (61) becomes, if h is constant as usual, 

/** p» /m 

p =s I I \ {,a cos x-\- t bsmx cos j/ + fi sin a; sin j/)u'sin x.dx dy du 

-5- I I I w 2 sina;.<fo<%c?u, 

in which equation the limits are so taken as to include a shell bounded by 
two similar ellipsoids, the plane t a]b, a plane through t a inclined to the 
plane of t a and t b at an angle tan -1 [(c-hb) tan y\ and a cone with its vertex 
at the origin and a section of the ellipsoid by a plane || to the plane of t b 
and fi for a directrix, the distance of this last plane from the origin being 
,a cos x. By integration we find 

3 1 — m 4 ,ay sin 2 x-\- ,h sin y(x— sinrccosa;)-)-^! — cosy) (a; — sin a; cos a;) 

iP» — 4 • ;PZ^s ~ 2^/(1— cos x) 

For one-eighth part of the ellipsoidal shell x = y = \u, and 

For one-fourth x = n, y = %it, 

If we had integrated from — y to + y instead of from to y the result 
above would have been unchanged except that the term containing c would 
not appear; if then x be taken from to n we shall have a lune lying be- 
tween two planes through a making equal angles with the plane of a and 
jb, of which the c. g. will be given by 

_ 3?r6 s'my 1 — u* 
' p o~TQ"~y~^'T=u*' 
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To obtain the solids reaching to the center we have only to make u = 
in each of the above equations. For the surface of the ellipsoid we have 
1— m 4 = l+u+u 2 +u s _ 4 

1— «3 l_|_ u _|_ tt 2 3' 

where « = 1 ; and this value substituted in each of the above expressions 
will give the desired results. Of course we have only to make a = b — o 
to get the corresponding results for the sphere. 

As an example of the use of (64) let us apply it to the torus generated by 
the revolution of the circle ,p = bj + a(i cos t -\-j sin t) about the axis of i. 

We shall have then e = i and for the equation of the torus 
t p = i th + n [ai cos t -f- (b + a sin t)f]i~" ri ' 1 '. 

Ve<p = Vi [ai cos t -f- (b + « sin i)J] = (6 + a sin i) A, ^\t) = a( — i sin £ 
+y cos £), Tip' Ve<p = o(6+a sin <), $s y p = $1,^ = — a cos t. Therefore 

/ir— < f— 2th ~ 1 +<ft 

iadcost — i "7c(6 + asin<) (6+asin<)cZi 
+ < L. _l -« 

+20 f* \b-\- a sin £)cft, 

The limits are taken from — to -f- 9 and from £ to n — t so as to get a 
portion of the surface symmetrical about j. From this we find on integr'n 
__ [b\n— 2t)+4ab cos t+%a\n—2t~2 sin* c os QJjsin fl 
/,0 ° ~" dJbn—2bt+2aeost) " ' 

from which by giving proper values to and t we may find the c. g. of any 
portion of the surface symmetrical about the planes of i and j, and j and k. 
We will next apply eq. (67) to the case of a hyperboloid of revolution of 
one sheet formed by revolving the straight line / p= / ^(t)=:bj -\-ts about the 
axis of i. For convenience we will suppose Sej = 0, so that the generatrix 
in its initial position meets / at a distance b from the origin and is J_ taj. 
For the equation of the surface we have then ,p = i' h ' hlT (bj-\-ts)i~ th '^' r , i 
replacing the e of eq. (67) .-. (p\t) = s, Vijp = Vi(bj + ts) = blc -f tVie, 
V/p/f>'= V(bj+te)e=bje, Stf / <p'Vi i <p==bSje{bk+t,Vie)=-b 2 8ie, Sitf=tSis. 

... /j0o=s |.i=!Lj dti-iSis-i^ibk+tVie)^ dt-i-Ojj dt 

^■S[ Wl(W1 ^( ,J,J,_1 ) [M(W+1{WF " ] ] 

=|-^[*(« a +«ir , ^-*r 1 (*' 1Sl -i)P*+l(*»+<i)^]]- 

Suppose <! = 0, u = and X = tz, for a half revolution, 

••• <P* = %i- 1 Siet 2 +—Viet 2 + ^-k 



